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Relaxation of a 1-D gravitational system 
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We study the relaxation towards thermodynamical equilibrium of a 1-D gravitational system. 
This model shows a series of critical energies Ecn where new equilibria appear and we focus on 
the homogeneous (n = 0), one-peak (n = ±1) and two-peak (n = 2) states. Using numerical 
simulations we investigate the relaxation to the stable equilibrium n = ±1 of this A^— body system 
starting from initial conditions defined by equilibria n = and n — 2. We find that in a fashion 
similar to other long-range systems the relaxation involves a fast violent relaxation phase followed by 
a slow coUisional phase as the system goes through a series of quasi-stationary states. Moreover, in 
cases where this slow second stage leads to a dynamically unstable configuration (two peaks with a 
high mass ratio) it is followed by a new sequence "violent relaxation/slow coUisional relaxation" . We 
obtain an analytical estimate of the relaxation time t2^±i through the mean escape time of a particle 
from its potential well in a bistable system. We find that the diffusion and dissipation coefficients 
satisfy Einstein's relation and that the relaxation time scales as Ne^'"^ at low temperature, in 
agreement with numerical simulations. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The thermodynamics and dynamics of systems with 
long-range interactions have been the focus of many stud- 
ies in recent years [l|, Q- Indeed, such systems exhibit 
many peculiar features due to their long-range nature and 
to their non-additivity, such as ensemble inequivalence 
between micro-canonical, canonical and grand-canonical 
ensembles and regions of negative specific heat Q- On 
the other hand, their dynamics presents many interesting 
phenomena *4\. In particular, the relaxation to thermo- 
dynamical equilibrium can be very slow and diverge with 
the number N of particles [3, la, 13 ■ Moreover, this re- 
laxation often proceeds in two steps with very different 
time-scales. First, there is a coUisionless relaxation over a 
few dynamical times which involves collective dynamical 
instabilities (this step is also called violent relaxation in 
astrophysics (g)- Secondly, after the system has reached 
a mean-field equilibrium a coUisional relaxation associ- 
ated with two-body encounters, or more generally due to 
the discrete nature of the matter distribution which gives 
rise to fluctuations with respect to the smooth mean-field 
potential (finite N effects) , leads to a slow relaxation to- 
wards statistical equilibrium over a time-scale which di- 
verges with N Q. 

A prototype of such systems with long-range interac- 
tions is the Hamiltonian Mean Field (HMF) model de- 
fined by a cosine interaction for particles moving on a 
circle |3|- It has been shown [3 that this system first 
converges to a stable stationary solution of the mean- field 
Vlasov equation. Next, the relaxation to thermodynam- 
ical equilibrium proceeds over a much longer time-scale 
through a slowly varying sequence of stable stationary 
states of the Vlasov dynamics. This process implies that 
the dynamics of the system strongly depends on the ini- 
tial conditions since there are an infinite number of stable 
stationary solutions of the Vlasov equation j^, [i3 ■ 



In this article we study the relaxation of the 1-D grav- 
itational system described in details in [ll|. This One- 
dimensional Static Cosmology (OSC) model consists of 
particles moving between two reflecting walls within an 
external potential V which balances the 1-D gravita- 
tional self-interaction $ so that the homogeneous state 
(i.e. constant density) is an equilibrium solution. This 
model also corresponds to the evolution of 1-D density 
fluctuations in a 3-D cosmological background over time- 
scales much smaller than the Hubble time so that the 
expansion of the universe can be neglected. As shown 
in 11], this system exhibits a series of critical energies 
Ecu [Eel > Ec2 > •■)• At high energies the only stable 
thermodynamical equilibrium is the homogeneous state 
(called "n = 0" ) and at each transition Ecn two new ther- 
modynamical equilibria (±n) appear. Moreover, equilib- 
ria ±1 are stable below Ed (while the homogeneous state 
turns unstable) whereas other equilibria ±n with n > 2 
are unstable (both from a thermodynamical and a mean- 
field dynamical analysis) except for equilibrium ri = 2 
which becomes stable for the Vlasov dynamics at low en- 
ergy but remains thermodynamically unstable. In this 
article we study the relaxation of the OSC model start- 
ing from either the n = or n = 2 equilibria, which 
allows us to investigate both coUisionless and coUisional 
processes. We first recall in sect, mthe thermodynam- 
ical properties of the OSC model. Next, we study the 
relaxation of the system, starting from the homogeneous 
state in sect. Illl'XI and starting from equilibrium n = 2 
in sect. IIII 51 Finally we conclude in sect. If VI 



II. THE OSC MODEL 
A. Description of the model 

The OSC model [lj| consists of N particles of mass 
771 which move along the a;-axis in the interval < 
X < L (with reflecting walls) within an external concave 
quadratic potential V{x) and which interact through 1-D 
gravity. Thus, the Hamiltonian Tijv of the system is: 
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where Vi = it is the velocity of particle i (we note by 
a dot the derivative with respect to time t), g is the 
coupling constant of the 1-D gravitational interaction and 
the external potential V{x) is: 



V{x)^~g-p[{x-L/2'f + Ly4], 



(2) 



where p = M/L is the mean density of the system (M — 
Nm is the total mass). In [ij the thermodynamics and 
stability properties of the OSC system were studied in the 
mean-field limit (i.e. continuum limit) where the mass m 
goes to zero at fixed density p. Then, the gravitational 
self-interaction is described by the potential <i>(x) with: 



$(a;) = 5 / dx'p{x') \x - x'\, 



(3) 



and the dynamics of the system is governed by the Vlasov 
equation for the phase-space density f{x,v,t) within the 
total potential (jj: 



dt dx dx dv 



<^ + V. 



(4) 



The main feature of the OSC model which comes from 
its cosmological context is the appearance of the back- 
ground potential V{x) which ensures that the homoge- 
neous configuration p = p is a solution of the equations 
of motion (this is the counterpart of the Hubble flow for 
the expanding universe). Alternatively, the constant den- 
sity p defined from the potential V can be interpreted as 
a cosmological constant (in which case it is not neces- 
sarily equal to the mean matter density) if we work in 
physical coordinates. The statistical mechanics of such 
a system has been studied in the 3-D case in [l^, Il3| . 
Here, following [ll| we interpret the OSC model as de- 
rived from a simple cosmological framework without cos- 
mological constant. Then, the external potential V ap- 
pears through the change to comoving coordinates [iJI 
and the Hamiltonian |^ applies on time-scales which are 
shorter than the Hubble time (so that the expansion of 
the universe can be neglected and there is no explicit time 
dependence in the Hamiltonian (^). However, for our 
purposes the OSC model can also be studied for its own 
sake as a simple model of systems with scale-free long- 
range interaction which exhibits an interesting behavior 



(in particular the series of critical energies discussed in 
[ll| and recalled below in Ec ljlll) '). In the following we 
shall decompose the total energy E from ^ into its ki- 
netic (K), self-gravity ($) and potential (V) components 
as E = K + '^+V with: 
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B. Thermodynamical equilibria 

At thermodynamical equilibrium the phase-space den- 
sity is the usual MaxwcU-Boltzmann distribution: 



f{x,v) = p{x)d^e-^-"'^ ex e-''["V2+0(.)]^ 



(7) 



where we introduced the inverse temperature [3 — 1/T . 
The density is related to the potential by the Poisson 
equation modified by a constant term due to the back- 
ground V: 



-— ^ — 2g{p — p) with (j) = if + If and p — pe ^'^. 
dx'' 



(8) 



Here we introduced for convenience the offset Tp defined 
by the condition p = p for 1^9 — 0. Therefore, the ther- 
modynamical equilibrium is set by the equations: 
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1) and (p'{0) = ip'{L) = 0. (9) 



The homogeneous state p — fj (i.e. (p = 0) is a solution of 
Eqs.lO which we label as "equilibrium n = 0". Its total 
energy is: 



Ea = 



MT gM^L 



(10) 



Note that the energy Eq is bounded from below by 
£^min(0) = -gAPL/G. Besides, as shown in [ill, Eqs.® 
yield a series of critical temperatures Ten and energies 
Ecu [n = 1,2, ..) where new thermodynamical equilibria 
appear, with: 
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Thus, at high temperatures above Td the only thermo- 
dynamical equilibrium is the homogeneous solution p — JJ 
{n — 0). Below Td this state becomes unstable and two 
new stable equilibria "n = ±1" appear. The state n — 1 
corresponds to a density peak at a; = and a density 
minimum at x = L, whereas the state 77 = —1 is its 
reflection through x = L/2. Similarly, at each critical 
temperature Ten two new equilibria ±77 appear, which 



consist of n half-oscillations (i.e. from a density peak 
to a density minimum). The state n > shows a peak 
at a; = whereas the state n < shows a minimum at 
X = 0. As described in 11] the OSC model can also be 
extended to the whole real line by symmetry with respect 
to a; = and periodicity of 2L. Then, the state n < is 
a mere translation of the state \n\ by L/\n\. Thus, the 
equilibrium n — 2 displays a density peak at each bound- 
ary X — 0, L and a minimum at x = L/2 (whereas state 
n = —2 shows one density peak at x = L/2 and two 
density minima at x — 0,L). Moreover, all thermody- 
namical equilibria ±n can be obtained from state n = 1 
at a rescaled temperature: 



Pn(C;/3) = Pi(C;4)' P-n(C;/3) = Pi(C--;4)' (12) 



where we wrote explicitly the dependence on inverse tem- 
perature (3 and we defined the dimensionless coordinate 
Cby: 



X ^ X . , ^ 1 , L _„, 



where we introduced the ratio Cl of the system size L 
to the Jeans length Lj. All equilibria ±n with n > 2 
are thermodynamically unstable (in the three micro- 
canonical, canonical and grand-canonical ensembles, see 
|llj|'). They also exhibit a dynamical linear instability in 
the continuum limit where the motion follows the Vlasov 
equation, except for state n = 2 which becomes linearly 
stable very close to Tc2- Therefore, above Td the sta- 
ble thcrmodynamical equilibrium of the discrete system 
(Q is the state n — with p — 'p whereas below Td it 
is the state n — 1. However, as discussed in quasi- 
stationary long-lived states can also exist and below Tc2 
we can expect the equilibrium n = 2 to be such a long- 
lived configuration. 

At low temperature (C^ -^ oo) the density profile of 
the peak in the equilibrium n = 1 obeys ^j : 
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(14) 



whereas the density minimum at x — L scales as: 

Pi(a)-pe-<^i/2 = pe--^^-/2^. (15) 

For other states ±n we can obtain the asymptotics from 
p2|l . Note that the minimum energy which can be 
achieved by the OSC system is ii^min(l) = ~gM'^L/2 
when all particles are at rest at the same boundary, 
either x = Q ov x — L. It can only be reached by 
the equilibrium n = ±1 at zero temperature whereas 
other equilibria ±n with n > 2 have energies above 
E^^{n = 2) = -gAPL/4. 



C. Numerical simulations 

From Eq. QJ we obtain the equations of motion of the 
N particles as: 

X, - 2g-px, = g-pL f ' ^ ' - 1 j , (16) 

where N^ (resp. N~) is the number of particles located 
to the right (resp. to the left) of particle i, that is with 
X > Xi. In practice, following 15. 16J we rank the par- 
ticles in increasing order of Xi and when two particles 
"collide" we exchange their velocities so that the order- 
ing remains valid at all times: Xi < x^+i. Therefore, the 
numbers of particles to the left and to the right of parti- 
cle i are constant: N[' = i — 1 and N^ = N — i. Thus 
between two collisions the dynamics of particle i is given 
by: 



Xi{t) = x""^ + (a;° - x""^) cosh 
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td 



yn 



(17) 



where we defined x'^'^ as the (unstable) equilibrium posi- 
tion of particle i and idyn the typical dynamical time (for 
the homogeneous configuration): 
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In Ea. ((T7|) x° and vf are the particle coordinate and ve- 
locity at time i° (which is taken as the time of the last 
collision). We use the event-driven scheme of [l^ to fol- 
low the dynamics of the system. We store the position 
and velocity {x°,v^) of all particles as well as the time i^ 
of their last collision (initially we set i^ =0). Then, we 
compute the collision time of each particle with its neigh- 
bors and we store the results in a heap structure, so that 
the next collision is at the root of this heap. Next, we 
advance to this first crossing by evolving the two collid- 
ing particles with Ea. H17|l . exchanging their velocities at 
collision and updating their last collision time i°. Finally, 
we compute the new three collision times associated with 
these two particles and their neighbors which we store in 
the heap. At the next step we take care of the next 
collision. Therefore, we advance from one collision to the 
next and at each crossing we only need to update the tra- 
jectories of the colliding particles and their next crossing 
times. The collision time between particles i and i + 1 
can be obtained analytically from Ea. (|17|l . Indeed, this 
yields for the distance between neighbors an expression 
of the form Xi+i — Xi = A + Be*/*'^^" -I- Ce"*/*'^*'" where 
A, B and C are constants. Then the condition x^+i — Xi 
leads to a quadratic equation for the variable y — e*/*''^- 
which is easily solved. Therefore, both the trajectories 
and the collision times are computed "exactly" from an- 
alytical formulae. Their accuracy only depends on the 
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FIG. 1: (Color online) The transition time to->±i from the 
unstable homogeneous state n = to the stable equilibria n — 
±1 as a function of total energy E, below the critical energy 
Ed. At each energy the crosses are the numerical results 
obtained for 200 realizations of the initial condition n = 0, for 
a system of A'' = 50 (left panel) and N = 1000 particles (right 
panel). The solid line is the mean transition time obtained 
by averaging over these numerical results. The squares which 
are slightly shifted to the right show the transition times to 
states n = ±1 for realizations which happen to first relax to 
state n = 2 over a few dynamical times. The dashed curve is 
the average obtained for these systems. 



numerical accuracy of the computer and does not involve 
a discretization procedure to compute integrals or differ- 
ential equations. 



III. RELAXATION TO THERMODYNAMICAL 
EQUILIBRIUM 



A. Homogeneous initial state n = 



We study in this section the dynamics of the discrete 
OSC model (^ starting from initial conditions defined by 
the homogeneous thcrmodynamical equilibrium n = 0, 
i.e. Eq.Q with p = p. At high energies E > Ed 
the equilibrium n = is stable (and there are no other 
thcrmodynamical equilibria) and we checked numerically 
that the system remains in this configuration. Therefore, 
we focus here on systems with total energies E below the 
critical energy Ed where the state n = becomes un- 
stable (both from the thcrmodynamical analysis and the 
mean-field Vlasov dynamics) and two new stable equilib- 
ria n = ±1 appear, characterized by a density peak at 
X = (n = 1) or at a; = L (ri = —1). States ±1 are sym- 
metric with respect to x = L/2 and are essentially the 
same configuration. Hence we study here the relaxation 
of the OSC system from the homogeneous state n = to 
the thcrmodynamical equilibria ±1. 



1. Transition times 

We first investigate the dependence on the energy E 
and on the number of particles N of the transition time 
io^±i from state n = to either state n = ±1. Here 
we simply define the transition as the first time where 
\K - Ki\ < \K - Ko\ and |1 - ¥i| < |1 - ¥o| (i.e. the 
kinetic and self-gravity energies are closer to the levels of 
states n = ±1 than those of state n = 0). We display 
in Fig. ^ the transition time to^±i ^.s a function of the 
total energy E, obtained for 200 realizations (crosses) of 
the initial condition n = at each energy. Left (resp. 
right) panel corresponds to systems of iV = 50 (resp. 
iV = 1000) particles. The solid line is the average of 
these numerical results. 

We can distinguish three regimes in Fig. ^ First, at 
low energies below — l.TlScil the transition proceeds over 
a few dynamical times (crosses and solid line). Moreover, 
the comparison of both panels shows that this time-scale 
does not display a significant dependence on the number 
of particles. Indeed, it is set by the mean-field dynamical 
instability growth rate associated with the fact that the 
homogeneous equilibrium n = is linearly unstable for 
the Vlasov dynamics. 

Secondly, at energies above — 1.7|i?ci|, closer to the 
critical energy Ed, we find that the transition times 
to^±i can be quite large and increase with the number 
of particles. We can also note that they exhibit a very 
broad dispersion for different random realizations. 

Thirdly, at low energies below Ec2 — — 2.2|£'ci| we note 
that for some realizations of the homogeneous initial state 
the system does not directly relax to the equilibria n = 
±1. It first exhibits a transition to the equilibrium n — 
2 (which appears below the second critical energy Ec2) 
over a few tdyn and remains trapped in this two-peak 
configuration over a long time-scale until it eventually 
evolves to a one-peak state n = ±1. This leads to the 
very long transition times labeled to^2^±i shown by the 
squares in Fig. ^ (and the dotted curve for their average). 

We shall discuss these three regimes in more details 
below. 



2. Dynamical relaxation to n = ±1 

We first consider a typical configuration of the first 
regime which exhibits a relaxation to equilibrium ±1 
driven by the collective dynamical instability. Thus we 
choose a particular realization with N = 1000 particles 
of the initial condition defined by the unstable equilib- 
rium n = at _E = — 2|i?ci|. We show in left panel of 
Fig- El the evolution of the various contributions (0-® 
to the total energy E. The fluctuating curves are the 
gravitational self-energy $, the kinetic energy K and the 
external potential energy V of the iV— body system (from 
top to bottom). The constant curve labeled E is the to- 
tal energy which is conserved (as verified in Fig. |2J|. The 
other constant curves show the mean-field energy levels 
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FIG. 2: (Color online) Left panel: The evolution with time t 
of the various contributions K, $ and V to the total energy 
_E of a system oi N = 1000 particles. We display the curves 
obtained for a particular realization of the initial condition 
defined by the unstable equilibrium n = at energy E — 
— 2|_Eci|. The constant curves are the mean- field energy levels 
of K (solid lines), $ (dashed lines) and V (dotted lines) for 
the equilibria n = and n = ±1. The system starts at levels 
n = at t = and undergoes a transition to levels n = 1 
at i ~ 7tdyn. Right panel: The evolution with time of the 
masses Ml (dashed lines) and Mr (solid lines) located in the 
left and right parts of the system {x < L/2 and x > L/2). 
The constant lines show the values Ml , Mr of equilibrium 
n = l. 



of K, $ and V for the equilibria n = and n = ±1. Thus, 
we can see that the system starts at levels n = and dis- 
plays a transition to levels n = ±1 over a few dynamical 
times tdyn- The fact that the transition time-scale is of 
order of a few idyn and does not depend on the number 
of particles, as was checked in Fig. Qfoi' E < — 1.7|i?ci|, 
is due to the collective character of the instability. In- 
deed, as shown in [111] the thermodynamical equilibrium 
n = becomes linearly unstable below the critical energy 
Ed for the mean-field Vlasov dynamics. The growth rate 
of this instability saturates to e*/*''='" at low energies [ll| 
which sets the time-scale for the transition to equilibrium 
n = ±1. The right panel in Fig. |21 shows the evolution 
of the masses Ml and M^ located in the left and right 
parts of the system. It also clearly shows the dynami- 
cal instability which amplifies the small initial random 
imbalance Ml — Mr to build a left-peak state after a 
few idyn- We can note that the relaxation is not fully 
complete as the system still exhibits some mean devia- 
tion and oscillations from the levels n = 1 for i^ and V 
(and for masses Ml and Af/j), although the largest en- 
ergy component $ relaxes to its equilibrium level after 
~ 60idyn- We found that these small departures persist 

up to lO'^idyn- 

We present in Fig. 13 two snapshots of the density dis- 
tribution for the realization used in Fig. |21 at times t — 8 
and 12idyn- We can see that the matter distribution has 
indeed evolved from the homogeneous state to a one-peak 
configuration close to equilibrium n = I at t — Sidyn- 
However, it has not yet fully relaxed and some collective 



FIG. 3; (Color online) Two snapshots of the density distri- 
bution p{x) at times t — 8 and 12idyn. The histogram shows 
the matter distribution of the A'^— body system over 30 bins. 
The dashed curves are equilibria n = (constant density) and 
n — 1 (peak at a; = 0). 
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FIG. 4; (Color online) Two snapshots of the velocity distribu- 
tion f(v) at times t — 8 and 12idyn, in units of u/udyn where 
we defined «dyn = L/tdyn- The histogram corresponds to 
the TV— body system whereas the dashed curves are equilibria 
n — (narrow Gaussian) and n = 1 (broad Gaussian). 



oscillations persist as already seen in Fig. |21 which also 
lead to oscillations of the density and width of the peak. 

We show in Fig. 0| two snapshots of the velocity distri- 
bution at the same times t — 8 and 12tdyn- It has again 
evolved from the distribution n = to the distribution 
71 = 1 at time t = Stdyn- Although there are some oscil- 
lations in the tails (as seen from the right panel at time 
12tdyn and the fluctuations of K in Fig.[2Jl we can see that 
the relaxation of the velocity distribution is rather effi- 
cient and always remains close to the equilibrium n = 1. 

We checked that for a small number of particles A^ = 50 
the behavior of the system is the same in this energy 
range. The states n = ±1 and n = are well separated 
and the system exhibits a transition to equilibrium ±1 
over a few dynamical times. 
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FIG. 5: (Color online) The evolution of the energy compo- 
nents K, $, V and E (left panel) and of the masses Ml, Mr 
(right panel) as in Fig. |5| We show the curves obtained for 
a particular realization with N = 50 particles of the initial 
homogeneous state at _B = — 1.3|i5ci|. The system keeps wan- 
dering over states n = 0, ±1 and does not settle into a stable 
equilibrium. 
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FIG. 6: (Color online) The evolution of the energy compo- 
nents K, $, V and E (left panel) and of the masses Ml, Mr 
(right panel) as in Fig. |5| We show the curves obtained for 
a particular realization with A'^ = 1000 particles of the initial 
homogeneous state at -E = — 1.3|-Eci| which only relaxes to 
equilibrium n = 1 over ~ SOOOfdyn- 



3. Close to the critical energy Ed 

We now investigate the evolution for higher energies 
E > — 1.7|£'ci|, closer to the critical energy Ed- We 
first display in Fig. [S] the evolution of the energy and 
mass components for a system of iV = 50 particles at 
energy E = — 1.3|£'ci|. The right panel shows that the 
system evolves over a time scale set by the dynamical 
time tdyn, in agreement with Fig.^ but it does not settle 
into a stable equilibrium. Indeed, it keeps fluctuating 
indefinitely from left-peak to right-peak configurations 
and wanders over states n = 0,±1; this can also be seen 
from the fluctuating energy levels in left panel of Fig. [S] 
(we checked numerically that this fluctuating behavior 
remains unchanged up to lO'^idyn at least). Therefore, in 



this regime the transition time to^±i which was shown in 
Fig.^only corresponds to the first time where the system 
gets close to a state characterized by a density peak close 
to a boundary, as equilibria ±1, but the fluctuations are 
too large to let the system settle down in such a stable 
configuration. Hence it keeps exploring various states 
with both left or right overdensities. 

Next, we show in Fig. the evolution of a system of 
N — 1000 particles with the same energy E = — LSli^cil. 
Because of the larger number of particles the system 
starts with a more balanced state {Ml — Mj^)/M ~ 
l/vN and the fluctuations with time of various quan- 
tities are smaller. Therefore, we find that the system 
again wanders for a long time over states with Ml ~ Mj^ 
and even exhibits several oscillations from left-peak to 
right peak configurations but it now manages to even- 
tually settle down into a stable equilibrium n = ±1 af- 
ter ~ SOOOtdyn- Indeed, as the system starts close to 
Ml — Mfi small fluctuations are initially sufficient to 
evolve from Ml > Mli to Ml < Mb, but once it has 
converged close to either one of equilibria ±1 the fluctu- 
ations are no longer sufficient to escape to the symmetric 
state. 



4. Below the second critical energy Ec2 

Finally, we noticed in Fig.QJthat at low energies below 
Ec2 there are cases where the system does not exhibit a 
direct transition to stable equilibria ±1. It first evolves 
to the two-peak equilibrium n — 2 over a few dynamical 
times and next relaxes to the stable equilibrium n = ±1 
over a much longer time-scale. We show in Fig. [3 four 
snapshots of the phase-space distribution /(x, v) for such 
a realization at -E = — 2.5|-Eci|, with N — 1000 particles, 
which first relaxes to state n — 2. Thus, we can see 
that the system has already evolved from an homoge- 
neous configuration to a two-peak state at i = 6idyn- It 
remains in such a configuration until ~ 2250idyn while 
particles slowly diffuse out of the right density peak and 
at i = 2300tdyn it has relaxed to the stable equilibrium 
n — 1. Of course the energy levels show the same tran- 
sitions from n = to n = 2 and finally to n = 1. In 
agreement with the calorific curve shown in Fig. 4 of [I^l 
we note that for a fixed total energy E the temperature 
grows (the velocity distribution becomes broader) as we 
go from state n = to n = 2 and finally to n — 1. 
We shall investigate in more details the transition from 
equilibrium n = 2 to state n = 1 in sect. IIII Bl below. 
For systems with A^ = 50 particles we found that 58 
out of 200 realizations exhibit this two-stage behavior at 
E = — 2.5|-Eci| whereas for N — 1000 particles this only 
occurs for 30 out of 200 realizations. Indeed, as seen 
in [ij] the linear growth rate of the two-peak instability 
obtained for the mean-field Vlasov dynamics is smaller 
than the growth rate associated with a one-peak pertur- 
bation (higher wavenumbers are less unstable thanks to 
the finite temperature, following the usual Jeans instabil- 
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FIG. 7: (Color online) Four snapshots of the phase-space dis- 
tribution f{x,v) at times t = 0,6,2200 and 2300idyn for a 
system of N = 1000 particles. For this particular realization 
of the initial homogeneous state n — the system first un- 
dergoes a transition to the equilibrium n = 2 over 6idyn and 
only relaxes to the one-peak state n = 1 at i ~ 2250idyn. 



ity). Therefore, systems with a larger number of particles 
which follow more closely the mean-field dynamics should 
more frequently evolve directly towards a one-peak state. 



B. Initial equilibrium n = 2 

We now study in this section the relaxation to the sta- 
ble equilibrium n = ±1 starting from the equilibrium 
n = 2. 



1. Transition times 

We first show in Fig. [S] the transition times i2^±i ob- 
tained for various realizations of the initial equilibrium 
n = 2, as a function of total energy, for systems with 
A^ = 50 (left panel) and N = 500 particles (right panel). 
We also plot the mean transition times derived from these 
realizations (solid line) and the theoretical prediction 
(dashed-line) of ea. (|63|l . As in sect. lIII A ll wc defined the 
transition as the first time where jiiT— ii'i| < \K—K2\ and 
1$ — $i| < 1$ — $2! (i-e. the kinetic and self-gravity en- 
ergies are closer to the levels of states n = ±1 than those 
of state n — 2). We can see that there is a large dis- 
persion from one realization to another for small systems 
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FIG. 8: (Color online) The transition time t2^±i from the 
two-peak equilibrium n = 2 to the one-peak stable equilibria 
n = ±1 as a function of total energy E. The crosses corre- 
spond to various realizations of the initial equilibrium n = 2 
at a given energy, for systems of A'^ = 50 (left panel) and 
A^ = 500 particles (right panel). The solid line is the mean 
transition time obtained from these realizations whereas the 
dashed-line is the theoretical prediction (16311 . 



{N — 50) but the mean transition time agrees reasonably 
well with Ea. H63|l . For larger N where discrete effects 
are less violent {N — 500) the dispersion is much smaller 
and we recover the theoretical prediction H63|l which is 
derived in the limit of large N (which allows a pertur- 
bative analysis). In particular, note the steep increase 
at low energies when the two narrow density peaks are 
separated by an almost void region and there is a very 
slow diffusion of particles out of the smallest peak until 
the system reaches the abrupt transition point discussed 
below in Fig. 1101 where a collective instability merges the 
smallest peak into the largest one. We further discuss 
the dependence on energy of i2^±i in Sect. IIII B3l below 
where we detail the theoretical calculation of the transi- 
tion time. 

Next, we present in Fig. |5| the mean transition time 
t2~t±i as a function of the number of particles iV, for 
various total energies. We compare the numerical results 
(squares) with the theoretical predictions (dashed lines) 
of Eq. H63I) . We can check that the agreement is reason- 
ably good and that the transition time scales linearly as 
N . In the limit N —^ oo the equilibrium n = 2 becomes 
stable, in agreement with the mean-field analysis of [Ij 
where it was shown that this state is linearly stable for 
the Vlasov dynamics. 



2. Diffusive relaxation to n — ±1 

We first show in Fig. ^| the evolution of the vari- 
ous contributions to the total energy, for a particular 
realization with N — 500 particles of the initial condi- 
tion defined by the unstable equilibrium n = 2 at en- 
ergy E = — 2.8|ii^ci|- We can see that the system dis- 
plays a transition from levels n = 2 to levels n = 1 at 
t ~ 4.7 X lO^tdyn- We can note that there is first a slow 
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FIG. 9: (Color online) The mean transition time t2-.±i as 
a function of the number A'^ of particles. The various lines 
correspond to total energies E = — 2.73|i5ci|, — 2.9|i5ci| and 
— 3.1|i5ci| from bottom to top. The squares are the numerical 
results (averaged over many realizations, solid line in Fig. (HJ 
whereas the dashed-lines are the theoretical prediction 1631 1. 
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FIG. 10: (Color online) Left panel: The evolution with time 
t of the various contributions K, $ and V to the total energy 
i5 of a system oi N = 500 particles. We display the curves 
obtained for a particular realization of the initial condition 
defined by the unstable equilibrium n = 2 at energy E — 
— 2.8|-Bci|. The constant curves are the mean-field energy 
levels oi K (solid lines), $ (dashed lines) and V (dotted lines) 
for the equilibria n — 1 and n = 2. The system undergoes 
a transition from levels n = 2 to levels n = 1 at f ~ 4.7 x 
lO^tdyn. Right panel: The evolution with time of the masses 
Ml (dashed lines) and Mr (solid lines) located in the left 
and right parts of the system {x < L/2 and x > L/2). The 
constant lines show the values Ml, Mr of equilibrium n = 1. 
The two linear solid lines starting from Ml = Mr = M/2 at 
t — are the theoretical estimate 16211 for the mass transfer 
between both density peaks. 




FIG. 11: (Color online) Four snapshots of the density distri- 
bution p{x) at times t = 0,46500,46800 and SOOOOtdyn- The 
histogram shows the matter distribution of the A''— body sys- 
tem over 30 bins. The dashed curves are equilibria n — 1 and 
n = 2. 



drift up to ^ 4.7 x lO^tdyn where there is a sudden transi- 
tion towards energy levels of n = 1. This is more clearly 
seen in the right panel which displays the evolution of 
the masses located to the left {Ml at a; < L/2) and to 
the right {Mn at a; > L/2) of the system. We can see 
that there is a slow transfer of matter from the right peak 
to the left peak until ^ 4.7 x lO'^tdyn- Then, there is a 
sudden jump where the mass ratio reaches its equilib- 
rium value after a few dynamical times. This suggests 
that there is first a slow evolution over a few thousand 
tdyn with a quasi-static exchange of matter between both 
density peaks which remain close to local thermodynam- 
ical equilibrium states for separate clouds until the ratio 
of both peak masses becomes so large that the system 
turns dynamically unstable and relaxes towards equilib- 
rium n = 1 over a few dynamical times idyn- Indeed, as 
shown in pTl | both statistical equilibria n = 1 and n = 2 
are stable at low temperature for the mean-field Vlasov 
dynamics (although only state n = 1 is thermodynami- 
cally stable). By continuity, this means that there is a 
threshold where the system leaves the basin of attraction 
of two-peaks configurations to enter the basin of attrac- 
tion of the one-peak equilibrium. Indeed, two-peak states 
close to the equilibrium n = 2 (i.e. with a mass ratio of 
both clouds close to unity) have negative stability eigen- 
values and remain stable with respect to the mean-field 
Vlasov dynamics. On the other hand, two-peak states 
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FIG. 12: (Color online) Four snapshots ol the velocity distri- 
bution f{v) at times t = 0,46500,46800 and SOOOOidyn- The 
histogram corresponds to the A''— body system whereas the 
dashed curves are equilibria n — 2 (narrow Gaussian) and 
n = 1 (broad Gaussian). 



with a high mass ratio can be seen as small perturba- 
tions of the one-peak equiUbrium n — I which is itself 
stable. Therefore, such two-peak states undergo a col- 
lective instability which leads to relaxation towards the 
stable equilibrium n — 1. For the system shown in Fig. 1101 
this transition between two different stability regions of 
the mean-field dynamics occurs at i ~ 4.7 x lO^idyn. 

The solid lines starting from Ml = Mn = AI/2 at 
t = are the theoretical estimate (162(1 which has been 
normalized by a factor 0.2 of order unity to match the 
slope at early times of the mass transfer. This was then 
used to predict the transition times shown in Figs. I8I9I 
for a broad range of total energies and particle numbers. 

We present in Fig. ^2 four snapshots of the density 
distribution p{x) for the A^— body system as compared 
with the mean- field equilibria n = 1 and n — 2. In 
agreement with Fig. 1101 we recover the sudden transition 
at ~ 47000tdyn from a two-peak state to a one-peak state, 
after the right peak has slowly lost about half its mass 
to the left peak. We can also check that the density 
distribution relaxes to the statistical equilibrium n = 1. 

In a similar fashion, we show in Fig. E| the velocity 
distribution f{v), which is also seen to relax from the 
equilibrium ri = 2 to state n = 1. We note that the 
relaxation of the velocity distribution seems to be more 
efficient than for the density distribution. Indeed, at i = 
46500tdyn it is still close to the Gaussian n = 2 whereas 



FIG. 13: (Color online) Four snapshots of the phase-space dis- 
tribution f{x,v) at times t = 0,46500,46800 and 80000tdyn. 



at t = 46800idyn it is already close to the Gaussian n = 1. 
We display in Fig.^]four snapshots of the phase-space 
distribution f{x,v). It shows again that the transition 
from n = 2 to n = 1 proceeds in two steps, with a slow 
diffusion followed by a sudden disruption of the small- 
est density peak. In particular, we can see in the upper 
right panel of Fig. El that the right peak still exists as a 
distinct object a,t t — 46500tdyn and has disappeared by 

t = 46800tdyn. 

Finally, we show in Fig.^lfour snapshots of the phase- 
space distribution around the transition to the one-peak 
state. We can see that what is left of the right peak at 
t ~ 46770idyn exhibits a strong deformation at 46774tdyn 
and is absorbed into the left peak at 46776idyn. This 
describes in details the second step of the transition as a 
collective dynamic instability which merges the remains 
of the smallest peak into the largest one over a few dy- 
namical times idyn- 



3. Estimate of the mean transition time t2—>±i 

We describe in this section how the mean transition 
times i2-*±i shown in Figs. ISO can be estimated by an- 
alytical means. 

The relaxation towards thermal equilibrium of 
A^— body systems such as the OSC model is often studied 
through the stochastic dynamics of a test particle in in- 
teraction with the rest of the system. In such cases, the 
test particle experiences both a systematic drift (such 
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FIG. 14: (Color online) Four snapshots of the phase-space 
distribution f{x,v) at times t = 46770,46774,46776 and 
46778idyii around the transition to the one-peak state. 



as the dynamical friction suffered by a high- vel ocity star 
moving through a cloud of low- velocity stars [la 13) and 
a diffusion, due to random encounters, which both grow 
linearly with time. Moreover, if the correlation times are 
small one can use a Markovian approximation so that the 
relaxation of the test particle in the thermal bath is de- 
scribed by a Fokker-Planck equation 17], which is fully 
parameterized by the friction and diffusion coefficients. 
One may then describe for instance the relaxation of the 
velocity distribution of the test particle and its approach 
to thermal equilibrium |2(]Ll2ll|. 

In the case we consider here, the system starts in an 
inhomogeneous mean-field equilibrium with two peaks at 
the boundaries. We model the slow loss of matter of the 
smallest peak, displayed in Fig.^| as a diffusive process 
which makes particles located in the smallest potential 
well escape to the other potential well. Since the escape 
time is shorter for the smaller potential well we neglect 
the flux from particles moving from the deepest well to 
the smallest one. This is consistent with the right panel 
of Fig. ^1 which shows that the lightest peak is almost 
steadily losing mass. Here we shall neglect the evolu- 
tion with time of the small density peak due to mass loss 
and we simply consider the bistable system n — 2 made 
of the symmetric double potential well (j>{x) which has 
two minima at the walls x = 0, L and a maximum at 
X = L/2. Then, we focus on the behavior of a test mass 
m orbiting in the left peak and we study the evolution of 
its energy e (whence of its orbit) due to encounters with 



other discrete particles. Computing the associated fric- 
tion and diffusion coefficients we obtain the average time 
T required to reach the barrier energy (j){L/2) — (j){0) and 
escape from the potential well through a Fokker-Planck 
dynamics. This allows us to estimate the transition time 
t2^±i ^ T to relax to the thcrmodynamically stable equi- 
librium n — 1. 

In order to separate the mean-field dynamics from the 
discrete effects which give rise to the diffusion of particles 
we write the Hamiltonian Q as: 

Hat = m(Ho + 7^/) (19) 

and we defined the mean-field Hamiltonian Tio by: 



N 



-± + <^o{x,) + V{x,) 



and the interaction Hamiltonian Tij by: 



ni=e'^ 



3>j' 



3 



<^0{Xj 



(20) 



(21) 



In Ea. (|2()|l the mean- field gravitational potential $o is 
given by Eq.(|3l where p{x') is the mean- field equilibrium 
distribution {Tj). In Ea. (|21|) we added a factor e'^* for 
the computation of perturbative eigenmodes and we shall 
consider the limit w ^- 0"^. Thus, Tig describes the mean- 
field dynamics whereas Tii describes the discrete effects 
which vanish in the limit N —^ oo. Therefore, we consider 
Tij as a perturbation to TYq and we apply a perturbative 
analysis to the dynamics of individual particles. Here 
it is convenient to work with the action-angle variables 
(J, w) defined from the Hamiltonian TCq which describe 
the motion of a particle of energy e = v'^ /2 + (j)o along 
its non-perturbed orbit in the equilibrium potential (j)o = 
^0 + V from position a;_ to x+ p^l2a|: 



J 



2tt 



da;^2(e- Mx)), 



(22) 



and: 



= f7 



dx 



^2{e-Ux')) 



with 



^.-. (23) 



Thus, TYq = T^o(Jj) only depends on the actions Jj. On 
the other hand, thanks to the periodicity of 2ti with re- 
spect to w of the non-perturbed orbits x{J, w) we can 
write the two-body interaction as: 



gm\x — x 



k.k' —~oo 



0fc,fc'(J,J')e*^'"'-"'"''\ (24) 



which defines the Fourier coefficients 0fe.fe'(J, J'). Be- 
sides, from Eq.Q the mean-field potential ^o{x) reads: 



Mx) = gJdJ'dw'fo{J')\x-x'\ (25) 

= — E /dJ7o(J')0M(^,'^')e'''",(26) 
m ^ — ^ J 



k— — oo ' 
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where we used the canonical change of variable dxdt; = 
dJdw and /o is the mean-field equilibrium distribution 
0. This yields: 



Hi = e" 



j^j' 



\Y.Y.^>'^^'^'^^jy 



i{kw — k'w') 



3:3' ^.k' 



--EE /dJ7o(J')'/'M(-/>'^')e^''' 



j k 



Then, the equations of motion become: 



(27) 



Jj=—^{no+Hi), w,^-^{Ho+ni). (28) 

We write the action-angle trajectories { J(i), w(i)} as the 
perturbative expansions J = j'"-* + J^^-* + J^'^-' + ... where 
J*^'^' is formally of order k over Hi. At zeroth-order we 
simply have: 

i(") = _^ ^ 0, wf^ = ^ = f7(jf ), (29) 



dw 







dJj 



which yields the mean-field equilibrium orbits: 
jj°^ = constant and wf ^ wf\fd) + n{jf^)t. (30) 



At first order we obtain: 



J. 



fW 



dHi 



dwi ' ^ 



(1) 



dn 
dj' 



r(l) 



dHi 



(31) 



where we can substitute the zeroth-order orbits in the 
r.h.s. Using the property (j>k',k{J', J) = (t>k,k'{J, J')* ob- 
tained from Ea. (|^ a simple calculation yields (e.g. [23): 



with: 



7(1) 



j¥i' 



^ w^'^ - ^ (32) 



dw 



dJi 



1 ■'^■' 



l)k,k' 



^i{kw — k'w') 



2 ^^ ^^ uj + i(kn - k'n') 

jd' k,k' 



27r 
m 



5:5:/dj'/o(j') 



4'kfi _^ikw 



] k 



Lo + ik^ 



(33) 



At second order we need to follow the first-order orbits 
in the r.h.s of Ea. (|28() and we have for the actions: 



jf ^ = - 



E 






(34) 



j' 



Using again (f)k',k = 4'-k,-k' — 4>k k' ^^'^ taking the av- 
erage (with the equilibrium distribution fo{J)) over the 



actions and angles {J^^\ ^(''•'(O)} of other particles yields 
at order 1/A^ (note that (f>k,k' oc m and m ex 1/N): 



(j(2)) = 



dJ 



27r 
m 



fdj'MJ')J2 

•' k,k' 



^k''\4>k,k'? 



to^ + {kn - k'n')^ 



47r^ /■ t^ — T ujk'^(bk n(t''i n 



Nm? 



•^ k,k' 



ujkk'\(f>k,k' 



j"^ + (m - k'^'Y 



W(35) 



Then, using \\ni^^Q+ uj/{uj^ + x^) — i:5d{x) where 5d is 
Dirac's distribution, the limit cj — > 0^ gives: 



(j(2)) 



27r2 
m 



|a//„,r,i;(4-.'A 



k,k' 
|2,.!r /7,o iJmi\ 



x\(j)kM'rkSDikn~k'n'), 



(36) 



which does not depend on time. From these results we 
can obtain the mean drift and diffusion of the action J of 
a test particle. First, the average change (A J) = ( J(t2) — 
J(ii)) of the action over a time- interval At = t2 — ti reads 
at order 1/A^: 



,AJ. ,AJ(i'+Aj(2) 



At 



At 



> = iJ^'^) 



(37) 



since (J*-^-*) is constant and (J*-^-*) — when we average 
[23 over the angles ^'■^•'(O) as seen from Ea. H32|l . Next, 
from Eas. H32|l - H33|l the mean-square change ((AJ)^) 
reads at order 1/A^: 



{{Ajf 



3^"'^- /dj'/o(j')E 

■' k,k' 



k^l't'k.k' 



uj^ + {kn - k'n')^ 



X (l + e 



2ujAt 



2e'^^*cos[(fcrj-fc'f7')At] 



2c.ti 47r2 /■,„,,„,,„,,, „„ V- '' <Pkfi9kfi 



^„,2 / dJ'dJ"/o(J')/o(J")E 
X (1 + e^"^* - 2e"'^* cos[fcf7At]) 
The limit a; — > 0^ now gives: 
AAJf 



^ Uj"^ + (fcf])2 



(38) 



At 



-— j dj' w)y^2e\(t>k,k 

•^ k,k' 



1 - cos[(A:^ - A:'rj')At] 
■ At{kn ~ k'Vl'Y 

-±^ I dJ'dJ"UJ')W')Y,'^k^kfi^ 
1 - cos[A:f]At] 

At(fcr2)2 ■ 



'fc,0 



(39) 



As seen from Eas. (|32ll - (|33f) the action J and angle w of 
each particle in a given realization (i.e. without perform- 
ing any averaging) are modified by discrete effects over 
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a time-scale which grows as v ^ and as expected in the 
Umit N —> oo we recover the mean-field dynamics. Since 
we consider the limit of large number of particles, which 
justifies the perturbative approach (|19|l . let us consider 
in Ea. ()39|l time-scales of order At ~ ^/N -^ oo. Using 
the equality limt„,oo(l — cos tx)/tx^ = irSDix) we obtain 
in this large- A'^ limit: 



the scalings: 



' At 



= — j dj' w)y^e\<i,k,k'?5o{k^^k'Q'). 

m / ^-^ 

■' k,k' 

(40) 
Thus, both the mean drift and diffusion of the action oc- 
cur through resonances between particle orbits (see also 
[23 )■ Moreover, they both grow linearly with time. We 
can note that [2J| also obtained a diffusion coefficient of 
the form 139|1 which decays as 1/i with an oscillatory 
behavior for the HMF model in the limit of low tem- 
peratures where a density peak appears. However, |24| 
approximated the core orbits by an harmonic oscillator 
with an unique frequency flo so that there was no integra- 
tion over Q' as in Eas. (|SS|l . (|^ through the dependence 
r2'(J'). In order to take the large-time or large-iV limit 
it is necessary to take into account the distribution of 
orbital frequencies fo{fl). Note that our approach also 
applies to any Hamiltonian system with two-body inter- 
actions, such as the HMF model. At this point, it is 
convenient to change variable from the action J to the 
energy e( J) = u^/2-|-(/)o (both are defined from the mean- 
field Hamiltonian Tia). Using the expansion: 



Ae = f}AJ+l^(AJ)2 + 



(41) 



we obtain at order 1/A^: 



,A.^ 
'At' 



47r2 

771 



|d.77.(.7') E (4-^'^ 

•^ k.k' = l ^ 



x\(i3k,k'\'^kn6Dikn~k'n'), 



(42) 



and: 



,(A6) 



Stt 



At 



dJ7o(j') V \ct>k,k'\'k'n^D{kn-k'n'). 



k,k' = l 



(43) 
In Eas. l42ll - 143|l we changed the sums over k and k' from 
] — 00, oo[ to [1, (X)[ which gave a factor 2. We can check 
from Ea. (|42|) that the transfer of energy Ae(J' ^ J) = 
— Ae( J -^ J') is anti-symmetric over {J, J') so that the 
mean change of energy over all particles (with the dis- 
tribution fo{J)) vanishes. This is related to the conser- 
vation of energy by the Hamiltonian dynamics. On the 
other hand (Ae)^ is symmetric. We now focus on the 
low temperature regime T — > 0, Cl — > 00, where the ini- 
tial state n — 2 consists of two narrow density peaks 
at the boundaries x = 0,L. Then, from the asymptotic 
behavior (|14|l and the rescaling H12|) we obtain for the 
properties of the core which contains most of the mass 



Pc y.2 ^c 



^-2 



Qr 



n. 



dyn 



(44) 

where Pc, Rc,^c, Vc and Jc are the typical density, orbital 
radius, frequency, velocity and action of core particles 
and we defined fidyn — l/^dyn, Vdyn = L/tdyn and Jdyn = 
LVdyn, from the dynamical time tdyn of Ea. H18|) . Let us 
consider the dissipation rate 7 = —{Ae/At) of a halo 
particle with an orbit of order i/4 which probes a finite 
part of the left potential well (i.e. not confined to the 
small core at a; < Re). Its typical frequency, velocity and 
action are r2dyn,Wdyn and Jdyn- Since most of the mass 
(except for an exponentially small fraction, see Eas. (|14|) - 
(|15|l ) is in the core the dissipation rate obtained from 
Ea. (|42|) is dominated by core particles J' . Therefore, we 
have k/ J ~ C,l and k' / J' ~ Cl hence expression (|42ll is 
governed by its second term. Next, using: 



d/o (^, d/o , 



(45) 



for the equilibrium distribution and integrating by 
parts we obtain: 



47r2/3 



772 



dJ'W] 



00 

E I 

fc,fc'=i 



Cfe.fc' 



|2;^2o2 



k^n^5D{kVL^k'VL'). 

(46) 
As could be expected, we find that the energy drift (Ae) 
corresponds to a damping term (7 is positive) of order 
1/A^. Thus, the high-energy halo particle loses its en- 
ergy to lower-energy core particles. This is similar to 
the dynamical friction suffered by a high-speed star mov- 
ing through an homogeneous background of lower ve- 
locity stars [iM Il9l|. Moreover, the diffusion coefficient 
D — ((Ae)^/At) obtained in Ea. (|^ and the dissipation 
coefficient 7 are related by: 



D 



27 



27T. 



(47) 



Thus we recover the usual Einstein relation (|47|l for a 
halo particle. Next, from the definition (|24|) we have: 

(pk^k' — giTL / 2 — \x — x'\cos{kw)cos{k''w'). (48) 

Jo ■"' 

If we approximate each orbit by its first harmonic x = 
Rsm{w/2) we obtain: 

0777 R 

(f>k,k' ^ --Z ^4m for A: > 1, A:' > 1, R' « R. (49) 

Taking into account the exact orbital trajectory would 
remove the Kronecker factor 5k', 1 (higher k' would con- 
tribute) but would not change the scaling gmR''^ /R. 
Next, from Ea. H27|l we see that the action and angle of 
each core particle in a given realization (i.e. without per- 
forming any averaging) are modified by discrete effects 
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over a time-scale Ate of order: 



J 

Ate r 

J 



N 



whereas for halo particles we obtain: 
At/,, 



^dyni 



(50) 



^NCiUyn » At,, (51) 

where we used the scalings H44|l and (|49|) . Therefore, 
we see that the core particles are significantly perturbed 
over a time-scale which is much smaller than for the halo 
particles. Then, the resonances between halo and core 
particles will be detuned over times of order Ate during 
which the trajectory of the halo particle has only suffered 
small deviations. This means that the core particles act 
as an external noise characterized by small time-scales 
with respect to the halo particle. This justifies a Marko- 
vian approximation for the behavior of the halo particle 
and allows us to write a Fokker-Planck equation [T3 for 
the evolution of the probability distribution P{e) of its 
energy: 



dP 

'dt 



d_ 
de 



ilP) 



2 9e2 



(DP). 



(52) 



The orbit-averaged Fokker-Planck equation (|52|l clearly 
separates the slow changes in phase-space coordinates 
due to encounters from the fast orbital motion in the 
mean-field potential [l9|. Indeed, the dissipation and dif- 
fusion coefficients scale as 1/A^ from Eas. (|42|l - H43() . It is 
interesting to note that for the OSC model (^ particles 
which do not cross the test mass do not contribute. This 
can be checked from Eq. (|48l) by noting that in such cases 
the absolute value can be dropped which yields (f>k,k' = 
for fc 7^ and k' 7^ 0. This is related to the fact that 
in 1-D gravity the force is merely proportional to the 
relative number of particles to the left and to the right 
of the test particle, independently of distances, as seen 
from Eqs.Q or H16|l . Therefore, only particle crossings 
can lead to fluctuations of the force seen by the test mass. 
The stationary solution Pst of the Fokker-Planck equa- 
tion |S2) is: 



^st(e) « — e 



-(27/r>)e 



(53) 



We can note that this distribution agrees with the sta- 
tistical equilibrium state "n = 2" described by the 
Maxwellian (O thanks to Einstein's relation (|47|l . On the 
other hand, we must point out that the relaxation of core 
particles cannot be described by the Fokker-Planck equa- 
tion H52I) because there is no separation of time-scales as 
in Eas. H50|l - H51|l which prevents the use of the Markovian 
approximation. This also explains why the Einstein re- 
lation (|47|l does not hold for core particles (since we can- 
not any longer neglect the first term of Ea. (|42|l l. From 
Ea. (|52|l we can obtain the mean time r it takes for a 
particle starting at energy e_ in the potential well^n to 
reach the maximum potential energy e+ . It reads |17j : 



de e^(' 



de' 



Die') 



-Fie') 



(54) 



with: 



F(e) 



-m-'^' 



(55) 



Eas. H54|l - H55|) clearly show that the escape time is dom- 
inated by the time spent at large values of e, that is at 
large radii of order L/2 far from the small core. This 
yields: 



P^D 



,/3(e+-e_) _ ^ _ ^(^^ 



(56) 



On the other hand, the diffusion coefficient (|43|l can be 
simplified by first integrating over J': 



D=^y 



k.y 



dJ' 



dl7' 



/o(J')l'/'fc,fcffc'^'' 



(57) 



n'=kn/k' 



Since fl/ft' ~ 1/Cl "^ 1 the sum over k can now be 
approximated by an integral over J' which yields: 



D 



m ^- — ^ 



fc' = l 



dJ' fo(.J'Mk,k' 



, k'^n'' 



n 



(58) 



k=k'n'/n 



From the scalings (|44|l and 149|l we obtain: 



NCltdyn 



(59) 



Taking for the energy gap Ae = e+ — e_ the potential 
barrier A<I> — Acf) between the bottom of the potential 
well at X = and the maximum at x = L/2 we obtain: 



'P- 



tdynNCt 



-1-ln 



f3Aip] 

Pin ax 



Pn 



(60) 



where we used Eq.®. Thus, at low temperatures we 
obtain from Eas. H14() - (|15|l the asymptotic behavior: 



UynNCte'^"/'' ~ tdy„A^ 



T 



,7r='Tc2/2T 



(61) 



where we used the rescaling ()12(l . From Ea. (|6()|l we can 
estimate the mass flux from one peak to the other one 
as: 



— = — with — 

dt 2tAf tdyn 



0.2 



NCI 



1-ln- 



'mill Pmin 

(62) 
Here the factors M/2 and N/2 express the fact that only 
half of the total mass is within each peak whereas the 
factor 0.2 of order unity has been chosen so as to match 
the slope of the mass transfer at early times displayed in 
right panel of Fig.^| which shows the evolution of a typ- 
ical system with a transition time t2^±i — 4.7 x lO^tdyn 
equal to the average value obtained from 200 numerical 
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simulations (solid line in right panel of Fig.|Sl). Then, in 
order to take into account the acceleration of the mass 
transfer close to the transition we simply write the mean 
transition time t2-*±i as: 



h 



+±1 — 



AM 
M/2 



tM with AAf = {ML{n = 1) - M/2)/4, 

(63) 

where MLin ~ 1) is the final mass located at x < L/2 
in the equilibrium state n — 1 whereas the factor 4 ac- 
counts for the late steepening of the mass transfer esti- 
mated from Fig. E| The transition time l|63|l is shown 
by the dashed-lines in Figs. I8I9I We can verify that it 
agrees with numerical simulations and recovers the steep 
increase obtained at low total energies, as the exponential 
of the inverse temperature (see Eq.jniJ). The departure 
at higher energies close to the transition Tc2 is expected 
since the approach described above was performed in the 
low-tcmpcrature limit Cl S> 1 where the hierarchy H51|l 
holds. Close to Tc2 one cannot distinguish halo and core 
particles and one cannot use a Fokker-Planck equation as 
(jSSJ- Besides, Ea. l|61|l shows that the relaxation time in- 
creases linearly with the number of particles iV, in agree- 
ment with the numerical results shown in Fig. |5| This 
is different from the relaxation of homogeneous states 
obtained in the HMF model [2l| where the relaxation 
time was seen to grow faster than N (numerical simu- 
lations gave ~ A''-'^-^ [3, |23)- However, we must note 
that Ea. (|61|l was derived in the low-energy limit where 
the system is strongly inhomogeneous. Therefore, the 
physical process involved in the relaxation (the escape 
of particles from the smallest potential well) is rather 
different. On the other hand, [2g also found a relax- 
ation time proportional to N for the 1-D gravitational 
system (which is identical to the OSC model without 
the reflecting walls and the external potential V). How- 
ever, using the Fokker-Planck equation derived in |23 
for the diffusion in (u, a) space (velocity v and accelera- 
tion a) from a Markovian approximation, '26l obtained a 
Fokker-Planck equation in energy space as in Eq. H52I) but 
without the friction term. Indeed, as seen above the lat- 
ter arises from the correlations between particles which 
accumulate over time but these were not considered in 
|27| which computed the transport coefficients in (d, a) 
space in the limit of infinitesimal time-steps and next as- 
sumed a Markovian evolution. This procedure cannot be 
applied here since correlations do not decrease exponen- 
tially with time (see for instance the cosine dependence 
in Eas. lfS^ . lfS^ l and the dynamics cannot be described 
as Markovian over time-steps of the order of a few orbital 
times (which amounts to erase all correlations every few 
orbital times). Moreover, the dissipation term in Ea. H52() 
is clearly required by physical consistency to recover the 
stationary distribution (|53|1 and Einstein's relation. 

We can note that halo and core particles are rather 
strongly coupled in the sense that the dissipation and dif- 
fusion coefficients only decrease as a power of the ratio 
of orbital frequencies Vt' /Q, (through the factors {R' /R)^ 
in Eq.lgHJ and {il'/i})'^ in Ea. (|S5|l 'l which gives rise to 



the inverse power of (l in Ea. H59|l . This is significantly 
different from the coupling in 3D gravity between a small 
dense core and unbounded scattering particles. In that 
case the smooth gravitational interaction (for non-zero 
impact parameter) leads to a coupling and a dissipation 
which decrease exponentially with the ratio of frequen- 
cies 28] . In the case we consider here the large transition 
time observed at low temperature is not due to exponen- 
tially small couplings between halo and core particles. 
It is merely due to the usual Arrhenius factor e^^'^l^ 
associated with the diffusion through a ffiiitc potential 
barrier. Finally, we can note that our results also apply 
to any Hamiltonian system with a two-body interaction 
which shows similar equilibria states. We only need to 
use the relevant scalings over the temperature T (in our 
case written in terms of <^l). 



IV. CONCLUSION 

We have studied in this article the relaxation of a 1-D 
gravitational system (OSC model) which was originally 
derived from the formation of large-scale structures in 
cosmology Tl', 16, 2ff|. We have checked that the ho- 
mogeneous equilibrium state n = becomes unstable 
below the critical energy Er_\ and exhibits a relaxation to 
the one-peak state n = ±1 over a few dynamical times. 
This is consistent with the fact that the homogeneous 
state is unstable both from a thermodynamical analy- 
sis and from a dynamical mean-field analysis. Therefore, 
the linear instability of the Vlasov dynamics leads to a 
violent relaxation to the stable equilibrium n = ±1 (or 
to a nearby one-peak state) which develops through a 
collective dynamical instability. 

On the other hand, close to the transition Ed we have 
found that for moderate numbers of particles [N = 50) 
the fluctuations due to finite N effects are large enough 
to prevent the system from converging towards a stable 
equilibrium. Indeed, the system keeps wandering over 
left and right peak states since it can easily jump from 
one equilibrium to another one. For larger numbers of 
particles {N ~ 1000) the fluctuations due to discrete ef- 
fects are much smaller than the distance (in terms of 
energy levels or mass ratios) between homogeneous, left 
and right peak states (unless we go closer to Ed). Then, 
starting from the homogeneous unstable equilibrium the 
system wanders again for a long time (~ 5000idyn) over 
left and right peak states but it eventually manages to 
settle in a stable left or right peak equilibrium. 

Finally, at low energies below Ec2 we have noticed that 
in some cases the homogeneous state does not relax di- 
rectly to a left or right peak configuration but first con- 
verges over a few dynamical times towards a two-peak 
state close to equilibrium n = 2. Then, the system un- 
dergoes a slow coUisional relaxation towards a one peak 
state n ~ ±1. 

Next, we have investigated the relaxation to thermody- 
namical equilibrium of the two-peak equilibrium n = 2. 
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As expected, we have found that since this is a stable 
equiUbriuni of the mean-field Vlasov dynamics the relax- 
ation involves a slow diffusion process over a time-scale 
t2^±i which diverges with N. Moreover, after the small- 
est density peak has lost most of its mass the two-peak 
configuration (with a high mass ratio) becomes dynam- 
ically unstable and the system converges to a one-peak 
state over a few dynamical times. We have estimated 
analytically the mean transition time t2^±i by describ- 
ing the slow diffusion of particle energies due to finite N 
effects with a Fokkcr-Planck equation. We have found 
that the friction and diffusion coefRcients of halo parti- 
cles satisfy Einstein's relation and we have obtained a 
mean transition time of the form t2^±i ^ Ne^''^ which 
is proportional to the number of particles and grows at 
low temperatures. We have checked that this prediction 
agrees reasonably well with our numerical simulations. 
The relaxation involves an cfhcicnt coupling between the 
halo particles which extend close to the barrier at L/2 
and the core particles buried in the density peaks thanks 
to efficient resonances at high harmonics. Thus, although 
halo particles have a much smaller orbital frequency than 
core particles the dissipation and diffusion due to the 
build-up of correlations and encounters do not vanish ex- 
ponentially with the ratio of orbital frequencies. 



Therefore, we have found that the relaxation of the 
OSC model proceeds in a fashion similar to some other 
long-range systems. It first involves a violent relax- 
ation phase governed by dynamical instabilities, where 
the system converges to stable solutions of the mean- 
field Vlasov dynamics over a few dynamical times. This 
is followed by a second much slower coUisional relaxation 
phase where the system goes through a series of quasi- 
stationary states of the Vlasov dynamics until it reaches 
thermal equilibrium. Moreover, we have found that this 
slow evolution can be followed by another violent relax- 
ation step as the series of quasi-stationary states may 
lead to an unstable configuration (the least stable eigen- 
value increases along the series and eventually becomes 
positive) which quickly relaxes to a new quasi-stationary 
solution. Here we note that the slow relaxation of the 
system as it goes through the series of quasi-stationary 
states is due to dynamical constraints and not to metasta- 
bility as for the cases discussed for instance in 30] where 
the system is trapped in local entropy maxima. On the 
other hand, the route to thermal equilibrium clearly de- 
pends on the initial conditions. Therefore, the relaxation 
time exhibits a strong dependence on the initial condi- 
tions, as shown by the comparison of io-*±i and t2^±i 
studied in this article. 
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